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ANALYSIS OF HEAT TRANSFER IN A POROUS COOLED WALL 
WITH VARIABLE PRESSURE AND TEMPERATURE 
ALONG THE COOLANT EXIT BOUNDARY 
by Robert Siegel and Marvin E. Goldstein 
Lewis Research Center 

SUMMARY 

A general analytical solution is obtained for the heat -transfer behavior of a porous 
cooled wall of constant thickness. The fluid flowing through the porous wall is supplied 
from a reservoir at constant pressure and temperature. The boundary through which 
the fluid exits has specified variations in pressure and temperature along it in one di- 
rection so that the flow and heat transfer are two-dimensional. The flow within the por- 
ous material is assumed to be governed by Darcy's law; that is, the local velocity is 
proportional to the local pressure gradient. The local fluid and matrix temperatures 
are assumed to be equal so that a single energy equation governs the temperature distri- 
bution within the wall. The solution is obtained by transforming the energy equation into 
a potential plane where by virtue of Darcy's law the potential is proportional to the local 
fluid pressure. In this plane the energy equation is separable so that a general solution 
is obtained. Results are given for the normal velocity and the heat flux along the surface 
through which the fluid is exiting. Examples are carried out to illustrate the effect of 
variable pressure and the interaction of variable temperature and pressure along the sur- 
face. It is found that in many instances the behavior is close to being locally one- 
dimensional. 


INTRODUCTION 

The need to operate some portions of advanced power -producing equipment, engines, 
and research devices at high temperatures has provided interest in transpiration -cooled 
media. A porous metallic structure is formed by such means as sintering small metallic 
particles or rolling together and sintering layers of woven wire cloth. A coolant can then 


be forced through the porous matrix so that it exits through the boundary exposed to the 
high-temperature environment. When the porous boundary is being heated by a flowing 
stream, the transplant serves a dual purpose. It cools the metallic region as it flows 
through it; and the exiting coolant pushes away the hot fluid stream, thereby decreasing 
the heat transfer to the surface. Some possible applications are for cooling turbine 
blades, components in fusion reactors, arc electrodes, rocket nozzles, high-speed 
bearings, and vehicles reentering the Earth T s atmosphere. 

In many instances the coolant exits from the porous material into a region of var- 
iable pressure. Two examples of this are the pressure variations resulting from the 
large acceleration in the throat region of a rocket nozzle and the pressure variation 
around the surface of a turbine blade. Another example of this occurs when transpira- 
tion cooling is used to cool the nose of a body in high-velocity flow. The flow through a 
porous hemispherical shell in such an environment has been discussed in reference 1; 
the high pressure at the stagnation point produces coolant starvation in that region. 

When the pressure along the exit boundary is variable, the solution for the heat- 
transfer behavior in the porous medium becomes difficult even for a simple one- 
dimensional geometry. This is because the flow in the porous medium is two- or three- 
dimensional and can be quite complicated if the pressure is strongly varying. Since the 
velocity appears in the energy equation, the solution for the temperature distribution in 
the porous material then becomes difficult. The resulting equations can be solved nu- 
merically; reference 2 gives such a solution for a reentry vehicle nosetip in which both 
the thickness of the porous material and the external pressure are variable. There is 
little other literature available for these types of situations. 

In the present report an analytical method is developed for determining the porous 
cooling behavior for a wall of constant thickness with the exit pressure varying along it 
in one direction. An additional complication is also accounted for which will allow 
coupling the porous-wall heat transfer with the heat transfer of the external flow. This 
is that the surface temperature can have a prescribed variation along the surface in the 
same direction as the pressure variation. The final results relate the heat flux at the 
high-temperature surface to the imposed surface pressure and surface temperature 
variations. The relation between surface temperature and heat flux can be used as a 
thermal boundary condition for the energy equation governing the external flow. 

The analysis utilizes the fact that when the fluid within the pores of the material is 
in very good thermal communication with the surrounding solid matrix (and this situation 
occurs in many applications), then locally within the material the fluid and matrix tem- 
peratures are essentially equal. The heat transfer is then governed by a single energy 
equation composed of two terms, one representing the energy carried by the flowing 
coolant and the other the heat conduction within the matrix material. 

For the slow viscous flow that occurs in the pores for many porous cooling applica- 
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tions the velocity is governed by Darcy’s law, which states that the local velocity is pro- 
portional to the local pressure gradient. Since the velocity also satisfies the continuity 
equation, the pressure in dimensionless form can be regarded as the velocity potential . 
The porous wall can then be transformed into a potential plane. Since the coolant reser- 
voir is at constant pressure, the adjacent boundary of the medium is a constant potential 
boundary , while the potential is variable along the coolant exit boundary corresponding 
to the specified external pressure variation. The reason for considering the porous ma- 
terial in terms of a potential is that by transforming the energy equation into potential 
plane coordinates the general solution of the equation can be found by using separation of 
variables. By applying a special technique to adapt this solution to the variable-pressure 
boundary condition a general method for obtaining solutions to problems of this type is 
obtained. 


ANALYSIS 

Governing Equations 

Consider the porous wall of thickness 6 with constant effective thermal conduc- 
tivity k m (based on the entire cross-sectional area rather than on the area of the solid 
alone) and permeability k shown in figure 1. There is a fluid with constant density p, 
constant heat capacity Cp, and constant viscosity p flowing through the wall . Assume 
that the thermal conductivity of the fluid is very small compared with k m and that the 
pore size is so small that Darcy's law holds. Let u denote the Darcy velocity (local 
volume flow divided by entire cross section rather than by open area) of the fluid. We 
suppose that no changes occur in the direction perpendicular to the x-y plane so the 
situation is two-dimensional. 


y 



Figure 1. - Porous cooled wall. 
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Below the porous wall (see fig. 1) there is a reservoir which is maintained at con- 
stant pressure and temperature Pq and t^ , respectively. The pressure of the fluid 
along the upper surface of the wall p can vary in an arbitrary fashion along the wall as 
long as the inequality pg > p g is satisfied. Then the fluid flows from the reservoir 
through the wall and out through the top surface. Since Pg is constant, the fluid velocity 
at the lower surface will be in the y-direction. However, the direction of the fluid veloc- 
ity at the upper surface is not necessarily perpendicular to this surface. 

If the thermal communication between the fluid and the matrix is sufficiently good, 
the local fluid temperature will be approximately equal to the local solid-matrix temper- 
ature. We denote this common temperature by t. When these assumptions are made, 
the heat and mass flow within the porous material are governed by the equations given in 
reference 3. Thus 


C+* 

1 

o 

11 

O 

(1) 

u = -— Vp 

(2) 

P 


V • u = 0 

(3) 


Boundary Conditions 

As the fluid in the reservoir at t^ approaches the lower boundary Sg of the porous 
wall (see fig. 1), the temperature rises to tg, which varies in an a priori unknown 
fashion along the x-axis. In most instances of practical interest the fluid velocity is suf- 
ficiently high, and the fluid thermal conductivity sufficiently low relative to the thermal 
conductivity of the solid, that the temperature change from t^ to tg occurs in a thin 
fluid region near s Q . This region is therefore assumed to be locally one-dimensional; 
and since the velocity is perpendicular to the surface (as a result of the pressure being 
constant along Sg), there is no flow along this thin thermal layer. By applying an energy 
balance across this thermal layer the boundary conditions at y = 0 are found to be 


k 

m 


at 

9y 


ptyt - ^ > I u i 

► 


for y = 0 


p = Pq = constant 


(4) 
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On the upper surface s at y = 6, we shall suppose that both the temperature and 
pressure vary in some arbitrarily prescribed manner except that their asymptotic? values 
at x = -°° are each the same as those at x = -f*> . We shall denote this common asymp- 
totic value of the temperature by tj and the common asymptotic value of the pressure by 
Pj . Let tg - tj and pg - p^ be the maximum absolute values of the deviations of tem- 
perature and pressure away from and p^ along y = 6. Then the boundary conditions 
at y = 6 can be written as 

t = t x + (t 2 - t 1 )F(x/6) for y = 6 (5) 

P = Pj + (P 2 " P 1 )H(x/6) for y = 6 (6) 

where the functions F and H are less than or equal to 1 in absolute value and go to 
zero as x approaches ±°° . 


Dimensionless Equations 


The following dimensionless quantities are now introduced: 


X 2 


= p c p “(Pq - Pi) 


2k 


m 


N = 


V 1 ! 

H - too 


p = 


P2-P1 

Po’Pl 


(7) 

( 8 ) 

(9) 


5 


I 


X = - 


y = i 
6 


<p = 


P 


p 0 'Pi 


T = 


t ~ t cc 

t X -Ko 


u = & 


K (P Q “ Pi) 


U 


( 10 ) 


These dimensionless quantities are substituted into equations (1) to (3) along with the 
boundary conditions (4) to (6) to obtain 


V 2 T - 2AU • VT = 0 

(11a) 

S- 

II 

(lib) 

V • u = o 

(11c) 

— = 2X lulT for Y = 0 

0Y 

(12a) 

<p = 0 for Y = 0 

(12b) 

T = 1 + NF(X) for Y = 1 

(13a) 

(p = 1 - PH(X) for Y = 1 

(13b) 


where 
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The porous wall in the dimensionless physical plane is shown in figure 2. When equa- 
tion (lib) is used to eliminate U from the other two equations (11), we obtain 

V 2 cp = 0 (14) 


and 


V 2 T - 2\Vc p • VT = 0 (15) 

Since the surface Y = 0 is at constant pressure, the function cp is constant along this 
line and it follows from equation (lib) that 

|U| = ^£ on Y = 0 (16) 

SY 

Hence the boundary condition (12) becomes 

for Y = 0 (17a) 

3Y 3Y 

cp = 0 for Y = 0 (17b) 

Equations (14) and (15), which are related to the flow and energy equations, are to 
be solved subject to the boundary conditions (13) and (17). Note that since the flow is 
independent of temperature, the dimensionless pressure cp can be found independently 
of the temperature by solving equation (14) subject to the boundary conditions (13b) and 

Y 


<p -- (P(X>, T - T(X), ip = 0 S (X) 


1 1 3 

I 2 

ip * 0 -"'" 

Figure 2. - Dimensionless physical plane. 
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(17b). Once this has been done, we can substitute the result into equation (15) and solve 
this equation subject to the boundary conditions (13a) and 17(a). 


Solution of Flow Problem 

Equation (14) shows that cp satisfies Laplace's equation, and hence there must 
exist a harmonic function !// and an analytic function W of the complex variable 

Z=X + iY (18) 


such that 


W = i// + icp 


(19) 


By using the inversion theorem for Fourier transforms it can be verified that the 
analytic function W whose imaginary part satisfies the boundary conditions (13b) and 
(17b) is given to within an unimportant real constant by 


where 


W = 


Z 


r>z _cw.dk 

2iri J_ x sinh k 


C(k) = 



H(X)dX 


(20) 


(21) 


and the Cauchy principle value is to be used in evaluating the integral in equation (20) 
across the singularity resulting from sinh k at k = 0. 

Physically, the change in 4/ between any two points is proportional to the volume 
flow of fluid crossing any curve joining these two points. Hence if/ must vary between 
-°° and +°° as X varies between and +°° . The mapping given by equation (20) 
therefore maps the infinite strip in the physical plane which is occupied by the porous 
wall into a region in the potential plane (W -plane) such as that shown in figure 3 in the 
manner indicated by the corresponding numbers in figures 2 and 3. 
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Transformation of Boundary Value Problem for T into Potential Plane 

When the independent variables in equation (15) are changed from X and Y to the 
coordinates \j/ and <p of the potential plane, the equation becomes separable. It can 
be shown from the Cauchy -Riemann equations that (see ref. 4) the first term in equa- 
tion (15) transforms as 


/v 

V 


2 


T = 


j 9 2 T | 9 2 T j 

dW 

\9i// 2 dcp 2 / 

dZ 


2 


The second term is transformed by using the relation (see ref. 3) 


V<p • 


VT = 


9T 

dW 

dcp 

9Z 


2 


Upon substituting these equations into equation (15) we find that 


dip 2 dcp 2 d( f > 


( 22 ) 


Since cp is constant (and equal to zero) at Y = 0, the boundary condition (17a) trans- 
forms into 

— = 2\T for (p = 0 (23) 

dcp 


i 


9 


These equations can be transformed into a slightly more convenient form by intro- 
ducing the new dependent variable 


0 = e Mi-<p) T 


(24) 


Then equations (22) and (23) become 

+ sfe . x 2 e = o 

dll/ 2 d(p 2 


(25) 


— = X© at (f - 0, -oo < < +o° (26) 

0 (p 


Either by separation of variables (as in ref. 3) or by taking Fourier transforms 
(with respect to ip) the general solution to equation (25) which satisfies the boundary 
condition (26) is found to be 


9 = e M<P-l) + J_ 

2tt 



T r <P ~T a (p 

( T g + X)e ^ + ( Tg - X)e P 


(t^ + X)e P + (t p - X)e P 


dp 


(27) 


where we have put 


T 


0 




The first term on the right is the particular solution for uniform pressure and tempera- 
ture along the boundary Y = 1. Hence the integral term gives the deviation from this 
solution resulting from variable pressure and variable temperature along that boundary. 

The function D(/3) must now be determined so that the boundary condition (13a) is 
satisfied along the upper boundary 341 of the region in the potential plane. In order to 
introduce this boundary condition into the problem, put (fig. 2) 

iP s (X) ^ ip(X,l) (28) 
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It now follows from equation (20) that 


* S (X) = X + 



e ikX e- k C(k) & 
sinh k 


Since the real part of a complex number Z is (Z + Z*)/2, this becomes 





IkX e~ k C(k) 
sinh k 



e -ikX 

sinh k 


Upon changing the variable of integration from k to -k in the last integral and noting 
that equation (21) implies that C*(-k) = C(k), we get 


^ (X) = X + J- 
iiri 



e ikX e~~ k C(k) 
sinh k 



e ikX e k C(k) dk 
sinh k 


Hence 


^ S (X) 



-COsL^dk 

tanh k 


(29) 


Now \p is the value of ip along the boundary 341 in the potential plane and is a mono- 
tonically increasing function of X. Hence equation (29) can be solved for X as a func- 
tion of if/ g to obtain 


X = X(ip s ) 

It is convenient to introduce the functions h and f of if/ by 

s 

h (p s ) =H[% g )] 


(30) 


(31) 
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and 



( 32 ) 


Then equations (13a) and (13b) show that the values of cp and T along the boundary 341 
in the potential plane are given in terms of i// g by 


<P = 1 - Ph(i// g ) 
T = 1 + Nf(^ g ) 


for W on 341 


(33) 


It therefore follows from equation (24) that, on this boundary, 0 is given by 


XPh(^_) i— —I ^ — v 

0 = e 1 1 + Nf (^ s ) for w on 341 (34) 

Hence evaluating equation (27) for W on 341 and substituting in equations (33) and (34) 
gives along the boundary 341 

2 sinh XPh(i// ) + Ne s f(^_) 

S o 



Equation (35) is a Fredholm integral equation of the first kind which can be solved 
for D(/3) once the temperature and pressure distributions F and H (and therefore also 
f and h) have been specified. However, for numerical purposes it is much more con- 
venient to work with an integral equation of the second kind. In order to obtain an equa- 
tion of this type we add and subtract the term 


277 




d/3 
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in equation (35) to obtain 


XPh(^) 

2 sinh XPh (&) + Ne s 

s 



d/3 



r 


<- 


t r [ _Ph (^s) T ft 1 / \ " T fif Ph ^ T fi 1 

(t p + X)e ^ le ^ - lj + (t^ - X)e p | e ^ - lj 


(T p + X)e ^ + ( Tj3 - X)e ^ 


d/3 


Upon taking the Fourier transform of both sides of this equation and interchanging the 
order of integration, we get 


r(a) =D (a) + — 

27 T 



K(a-,0) D(/3)d/3 


(36) 


where we have put 


T(a) s T 


a 


r e-‘^[: 

, 7-00 


2 sinh XPh(i// ) + Ne 
s 


XPh(i//_ ) 


) 


W d^, 


(37) 


t ( T fi + X ) e ^ I_(«,/3) + (to - X)e ^1 (a, £) 

K(a;/3) = -2. — B B ± 


(38) 


( T/3 + *)e T P + ( Tj3 - X)e" T £ 


and 


V «,*. f e ‘ <e -“ ) ^[e iPT ^ > . 1 ] 

tA- 00 


dtf' 


(39) 


Although for most functions h and f it will be necessary to evaluate the integrals (37) 
and (39) numerically, these integrals are simply Fourier transforms and therefore can 
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be obtained almost instantaneously by using a "fast Fourier transform" computer 
routine. Equation (36) can then be solved numerically for D(/3) either directly or by 
iteration. 

The terms appearing in equation (36) are complex. It is sometimes convenient for 
numerical purposes to replace this equation by a real equation. To this end put 

E(a) =lliD(a) + i±l D *(a) (40) 

2 2 

Then E*(a) = E(aO and this shows that E(a) is real. On the other hand, it follows 
from equation (27) that 0 will be real only if D*(a) = D(-a) (see appendix A for veri- 
fication) . Hence 

E (a) = - - - D(o?) + - + - D(-o?) 

2 2 

and 

E(-a) =i±lD(a) +l-±D(-o l ) 

2 2 

Therefore 

D(a) = l±iE(of) + ^E(-a) (41) 

2 2 


Substituting this into the integral in equation (36) gives after changing the integration 
variable from /3 to -/3 in the second integral" 


T(a) = D(a) + 


/ fl±iK(a ;/ 3) +l^lK(a;-j8)]E(/3)d/3 

277 I L 2 2 J 

J _oO 


(42) 


Equation (40) can also be written as 


E(a) = ReD(a) + ^*D(a) 
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Hence, upon taking the real part of equation (42) and adding it to the imaginary part, we 
get after some manipulation 


ReT(a) + 9#ir(a) =E(a) + A. f jl[k(a;0)+K*(a;0)] -^[k(«; -0) -K* (a; -/3)]jE(/3)d/3 

t/-QO 


By further defining 


Y(a) = Re r(a) + 9m r(a) 

(43) 

k(o<; 0) = ReK (a; 0) + 9m K(cr, -0) 

(44) 


we get 


Y(a) 


= E(a) + — f k(ar;p)E(f3)d/3 
27 T J - oo 


(45) 


This is the desired real form of equation (36) . 


Computation of Surface Fluxes 


The physical quantities of principal interest in the problem are the normal velocity 

of the fluid leaving the upper surface v_ and the conduction heat flux crossing this sur- 

s 

face q . It is convenient to define the dimensionless normal surface velocity V and 
s s 

the dimensionless conduction heat flux Q by 

b 


v 5 v 

S « (P 0 - Pi) S 


(46) 


and 


Qc 




k m (t l - t ~) 


(47) 
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Then it follows from equations (10) and (lib) that 


v = 

s 3Y 


|Y=1 


( 48 ) 


Upon using the Cauchy-Riemann equations this becomes 


V 


s 


dxj/ 


0X 


Y=1 


and by using definition (28) 


V 


s 


d^ g (X) 

dX 


Substituting equation (29) now shows that 


V S (X) 



Je 2« e iW£dk 

tanh k 


The conduction heat flux into the surface q is given by 

s 


q 


s 




(49) 


By introducing the dimensionless quantities (10) and (47) this becomes 


Q = 21 
s 0Y 


Y=1 


Hence equation (24) shows that 


Q =e M<p-l) BO +XQ dy\ at Y = 1 


i0Y 0Y , 
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Substituting equation (27) for 30/3Y in this relation gives 
r 


3Y 


* e 

] + J_ l 

r 00 r 

1 

D(s) - 

i 

✓» oo 

J - / vp/ 

30 

T 

i 

1 

°(/ 3)/3 

(T fl + *)e 1 

2 ?ri 3 Y 1 

f ■ Jr - 

7 

X) 

( Tp + *)e 


T n<P 


~T 0 <P 


(Tp + X)e P + (j p - X)e P 


d /3 


- T ^ 


-T, 


eW d /3 at Y = 1 


Since the Cauchy -Riemann equations show that d\p/SY = -d cp/dX, it follows from equa- 
tions (13b), (24), (28), and (48) that this can be written as 


Q s e XPH ( X ) 



cosh[XPH(X)] 


+ Ne ^H(X) 



D(/ 3 )J_(X;/ 3 )e 


i/3* s (X) 




2ffi 



J (X;/3)e 
T /3 


W s (x) 


d/3 


(50) 


where we have put 


r [ 1 -PH(X)] -t [l-PH(X)] 

(t + X)e P ± (t - X)e p 

J ± (X;/ 3 ) = _fi P 

(t p + *)e ^ + (Tp - X)e P 


(51) 


Equation (50) can also be expressed in terms of the real quantity E(/3) by inserting equa- 
tion (41) and changing the variable of integration from /3 to -/3 in the integrals involving 
E(-/3) to obtain (note that J(-/3) = J(/3)) 
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Q ^ e XPH(X) = v s( x)( x | 2 cosh[XPH(X)] + Ne XPH(X) F(X)} 

+ — f E(|3)J (X;j3) 

2v 

/*°° 

- P — / E(J8) A J (X; /3)jcos 0 iWX) + sin 0 (X)ld0 (52) 

dX 277 Tp + L s S J 



ILLUSTRATIVE EXAMPLE FOR PRESSURE AND TEMPERATURE 
VARIATIONS HAVING GAUSSIAN FORM 

To illustrate the use of the analytical results some specific solutions will be ob- 
tained for some cases where the pressure and temperature vary in an increasing or de- 
creasing manner about a minimum or maximum value. This general type of pressure 
variation would be encountered, for example, when a slot jet impinges normally against 
a wall (ref. 5). In this instance the pressure at the surface is a maximum at the stag- 
nation point and drops off monotonically moving away from this point in either direction 
along the wall. In addition to its specific behavior, the results for the pressure and 
temperature variations that will be considered will yield some general findings that can 
be applied to other variations . 

The pressure variation along the surface is given by equation (6) and we let 


h/|^ = H(X) = e 




(53) 


where <r is a parameter that varies the rate at which the pressure changes from its 
maximum or minimum value pg at X^. For large X the pressure is p^. The mag- 
nitude of the pressure variation along the surface, relative to the pressure drop through 
the wall, is fixed by the parameter P. 
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The surface temperature distribution, given by equation (5), will have the same 
form if we put 


. F(X) . e- (X - X ‘ ,2 A* (. 

The temperature variation is thus symmetric about X^.; and its rate of variation is 
governed by o^, which can differ from the parameter in the pressure distribution. 
The range of temperature variation along the surface is specified relative to the tem- 
perature drop between the coolant reservoir and the hot wall by the parameter N. 

To obtain the normal exit velocity from the porous wall the function H(X) is in- 
serted into equation (21) . Upon performing the integration the function C(k) is found to 
be 


C(k) = -P>/jt 


-ikX -(kcr /2) 2 
cr e *e ^ 


(55) 


This function is then used in equation (49) to obtain the dimensionless velocity V„(X) 
along the surface. The integral can easily be evaluated numerically by using fast 
Fourier transforms. 

After V g (X) has been found it can be used in equation (52), which gives an expres- 
sion for the dimensionless heat flux being conducted into the porous wall. This equation 
also contains the specified pressure and temperature variations as given by P, H(X), 

N, and F(X). The function (X) appearing in the integrals is given by equation (29) in 

s 

terms of the known function C(k). The function E(/3) needed to evaluate the integrals in 
equation (52) is found by solving the integral equation (45) numerically (by using 
Simpson 1 s rule and then inverting the resulting matrix equations) in which the functions 
y(a) and the k(cv;/3) are given by equations (43) and (44). The functions T and K 
needed to evaluate equations (43) and (44) are given by equations (37) to (39) . The steps 
in the solution are indicated on the flow chart in figure 4. The integrals in equations (29), 
(37), and (39) can be easily evaluated numerically by using a fast Fourier transform sub- 
routine. In the solution of equation (45) it is necessary to replace the infinite limits on 
the integral by finite values. For the numerical results that will be presented, values 
of ±12 (or ±15 for a = 6) were found to be sufficiently large to yield results that did not 

r 

change when the limits were further increased. 
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Figure 4. - Equations used to compute surface heat flux. 


DISCUSSION 

Several illustrative sets of results are now discussed to demonstrate the influence 
of the parameters which appear in the pressure and temperature variations given in the 
previous section. In the first few figures there is a variable pressure along the porous 
surface through which the coolant is exiting, but the surface temperature is constant. 
Then in the remaining figures the simultaneous interaction of temperature and pressure 
variations is illustrated. 

One of the parameters arising in the analysis is X = pC p K(pQ - which 

gives an indication of the relative amounts of energy being transported by the coolant 
flow and by heat conduction. This can be demonstrated by considering a wall of uniform 
thickness 6 across which there is a uniform pressure drop Ap and a uniform temper- 
ature difference At. The energy gained per unit area by the coolant flowing through the 
wall is vpCp At and the heat flux conducted is -k m At 5. From Darcy's law 
v = -(k/V)( Ap/6) so that the ratio of convection to conduction becomes pC p tf Ap/k m |j. 
or 2A. In what follows, the pressure drop varies with position along the surface so that 
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X, which is based on p Q - Pp indicates the relative magnitudes of the parameter in the 
constant-pressure region at large |X | . To demonstrate the effect of various pressure 
and temperature variations the parameter X will usually be taken as unity. Then the 
effect of X will be examined for a variable pressure with the exit surface at uniform 
temperature and with the exit surface at variable temperature. 

The effect of a variable pressure of the form described by equation (53) is shown in 
figure 5, where the dimensionless imposed pressure difference is plotted in part (a) of 
the figure. For this example the maximum pressure difference p Q - p^ occurs at both 
large positive and large negative X so that the dimensionless pressure drop goes to 
unity in these regions. For the parameters chosen, the pressure difference decreases 
to reach one -half of its maximum value at X = 0. The parameter <x specifies how 
rapidly the pressure changes along the surface. 

The velocity leaving the surface in the direction normal to the surface is given in 
figure 5(b) . The velocity is nondimensionalized so that it is unity in the region of con- 
stant pressure at large |X | where the dimensionless pressure drop is also unity. Thus 
the normalization is such that, if the flow were everywhere locally one -dimensional, 
the dimensionless velocity would be exactly the same as the dimensionless pressure. 

The figure shows that this is essentially the case for the most gradual pressure varia- 
tion shown (Op = 6). For cr^ = 3, however, the normal velocity decreases at x = 0 to 
somewhat below the locally one -dimensional value of 1/2 as a result of the two- 
dimensional nature of the flow in the region of increased surface pressure. 

When the heat flux being transferred to the fluid exit surface is nondimensionalized 
as shown in figure 5(c), a locally one -dimensional behavior would produce values iden- 
tical to the pressure curves. This is evidenced by the fact that the dimensionless heat 
flux also becomes unity at large |x| . It is seen that even for the rapid pressure varia- 
tions considered, the heat transfer is very close to being locally one -dimensional . This 
figure shows the heat flux being transferred to the surface when the surface is main- 
tained at a constant temperature. Thus where there is low cooling (low exit velocity), 
the surface heat flux is low. Even though the normal velocities are a little below 1/2 
near X = 0, the heat fluxes do not go below 1/2; that is, the local heat flux does not 
diminish as much as the local normal velocity. This is the result of heat conduction, 
as will be demonstrated by figure 7. This figure shows that the heat flux does follow 
the velocity variation more closely when X is increased, thereby diminishing the rela- 
tive effect of heat conduction. 

Figure 5 showed the effect of the shape of the surface pressure variation on the exit 
velocity and heat flux along the upper surface. Figure 6 gives similar pressure, veloc- 
ity, and heat flux curves for a fixed shape of the pressure variation, but with various 
pressure amplitudes. Shown is the effect of both increased and decreased pressure drop 
near X = 0. There are only very small deviations from. locally one -dimensional be- 
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(a) Dimensionless pressure along exit surface. 



(b) Dimensionless normal velocity along exit surface. 
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x/6 

(c) Dimensionless heat flux along exit surface. 

gure 6. - Heat transfer and normal exit velocity along coolant exit surface for exit surface 
at uniform temperature t^. Op = 3; X = 1. 



havior even though the pressure changes from minimum to maximum values are rapid 
enough so that they occur only over a distance along the wall equal to about five wall 
thicknesses. 

The parameter X = pC^K( p^ - P^)/ 2k m /J is a measure of the relative magnitudes of 
cooling by the flow and the heat conduction within the medium. If we think of k m as 
fixed, then X changes in proportion to the pressure drop Pq - p^. Thus a small X 
provides decreased flow and a greater influence by the matrix heat conductivity. As a 
result of the smaller flow associated with a decreased X the heat flux at the surface 
would be expected to decrease. In figure 7(c) this effect on the level of the curves has 
been normalized out by dividing the dimensionless heat flux by X. Thus the shapes of 
the curves as a function of X show only the effect of the relative importance of the ma- 
trix conduction. An increased conduction (decreased X) permits a higher surface heat 
flux in the low-velocity region, as some of the energy can be redistributed by conduction 
within the wall interior. However, as shown by the curves, this effect is small for the 
present case where the exit surface temperature is constant. 

The remaining four figures illustrate the effect of simultaneous variations in both 
surface pressure and temperature. First, consider the special case shown in figure 8 
that is given by the curves for (pg - P])/(Pq - P-^) = 0. In this instance there is a temper- 
ature variation along the fluid exit surface as given by the uppermost curve in figure 8(a), 
but the exit pressure is constant. The temperature variation is given by the uppermost 
curve in the figure. The temperature difference across the wall (surface to reservoir) 
at large |X | is tj - t^, and the distribution shown increases to 1. 5 times this value at 
X = 0. Allowing the local surface temperature to increase will permit the local surface 
heat flux to be larger. Now examine the heat -flux curve for a constant coolant exit pres- 
sure (pg - Pj)/(Pq - Pj) = 0 in figure 8(b). It is found that this curve is practically the 
same as the temperature variation', that is, the heat-transfer behavior is almost locally 
one -dimensional . 

For all the cases considered up to now, an individual variation in either the surface 
pressure or the surface temperature has produced essentially a one -dimensional re- 
sponse in the surface heat flux. As a consequence, when both surface pressure and tem- 
perature are varied simultaneously, the heat flux should be approximately the product of 
the two imposed variations when all quantities are expressed in dimensionless form. 

This is borne out by the results in figure 8(b) for the nonzero values of (pg - p^)/ 

(Pq ~ Pj) • 1° this instance both the temperature and pressure variations are symmetric 
about X = 0. 

Figure 9 exhibits the same approximate behavior for a temperature variation that is 
not symmetric relative to the three pressure variations shown. Thus the region where 
the surface heat flux can be largest corresponds to the location where the product of 
pressure drop and surface temperature is the largest. 


24 



(a) Dimensionless pressure along exit surface. 




x/6 

(c) Dimensionless heat flux along exit surface. 


gure 7. - Heat transfer and normal exit velocity along coolant exit surface for exit sur- 
face at uniform temperature t 1 and various values of A. a p = 3; (p 2 - - Pj> = 1/2. 




(b) Dimensionless heat flux along exit surface. 

Figure 8. - Effect of magnitude of pressure variation on heat flux along exit surface for 
variable surface temperature. (t 2 - t 1 )/(t 1 - 1^) = 1/2; o p = 0| = 3; X p = X t * 0; A = 1. 
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(c) Dimensionless heat flux along exit surface. 


Figure 9. - Heat transfer and normal exit velocity along coolant exit surface for exit surface at variable tempera- 
ture t s (x) with three different pressure variations, (pj - Pj)/(Pq - Pj) ■ 1/2; (t 2 - tj)/(tj -t^i =1/2; o t = 3; 

X) = 2; X = 1. 


Figure 10 illustrates a situation where the temperature variation is displaced var- 
ious amounts from X = 0, which is the axis of symmetry for the pressure variation. 
Again as a reasonable engineering approximation the dimensionless heat flux responds 
as the product of the pressure and temperature variations. 

This response is altered somewhat as the parameter X is varied and this is shown 
in figure 11. If X is approximately 1 or larger, heat conduction has little effect rela- 
tive to convection and the results in figure 11 show a locally one -dimensional response; 
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(a) Dimensionless pressure and temperature along exit surface. 



(b) Dimensionless normal velocity along exit surface. 


r 
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(c) Dimensionless heat flux along exit surface. 


Figure 10. - Effect of various temperature distributions coupled with variable pressure on heat flux 
along coolant exit surface. (p 2 - pjl/tpg - Pj) = (t 2 - tjWttj - 1^) = 1/2; o p = o t = 3; X p = 0; X = 1. 




that is, for X = 1 the dimensionless q_(x) curve is essentially the product of the dimen- 
7 s 

sionless p (x) and t (x) variations. When the conduction is increased relative to the 
s s 

convection (reduced X), the heat flux is increased in the vicinity of the maximum tem- 
perature as a significant amount of heat can now be conducted away. 


CONCLUSIONS 

An analytical solution has been obtained for the heat -transfer behavior of a porous 
cooled wall of uniform thickness with variable pressure and temperature along the sur- 
face through which the coolant is exiting. The solution was obtained by transforming the 
energy equation into a potential plane, to obtain a separable equation. 

The use of the solution was illustrated by applying it to several combinations of var- 
ious surface pressure and temperature variations. This revealed a number of interesting 
features. In all cases the coolant exit velocity followed quite closely the imposed pres- 
sure drop across the porous wall . This was true even for situations where the pressure 
varied so rapidly that it went from its minimum to its maximum value in approximately 
five wall thicknesses along the surface. For more rapid pressure variations this locally 
one -dimensional behavior would eventually be violated, but such rapid variations would 
not usually be expected in practice. _ 

The parameter X is a measure of the ratio of heat convection to conduction. When 
X ^ 1 (relatively large convection) , the heat flux variation was found to essentially follow 
a variation given by the product of the imposed pressure drop and surface temperature. 
Hence the maximum q can be tolerated in the local region where both pressure drop and 
surface temperature are high. This locally one -dimensional behavior is quite important 
because it justifies the use of a locally one -dimensional analysis for engineering design 
applications . 

For X < 1, heat conduction begins to be of some importance and can alter the local 
one -dimensional dependence of heat flux on the surface temperature variation. For con- 
stant surface temperature the properly nondimensionalized surface heat flux changes 
slowly with X. With variable surface temperature there is a significant influence of X, 
and as conduction becomes more important (X decreases) the maximum heat flux shifts 
toward the location of the maximum surface temperature. 



The solution relating surface heat flux and surface temperature provides the bound- 
ary conditions necessary to couple the heat -transfer behavior of the porous medium with 
the convection in the external stream. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 27, 1971, 

132-15. 
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APPENDIX A 


CONDITION ON D(p) TO MAKE 0 REAL 


By putting 


T /3 ( P ~ T 0<P 

(t + X)e P + (t - X)e P 

-£ £ (Al) 

(t^ + X)e P + (t^ - X)e P 

equation (27) can be written as 

Q = e M<P~ 1) + J_ f D(/3)J(/3)e W dj3 (A2) 

2tt * / -°° 

o 

The function J(/3) is real; and because it depends only on /3 , it satisfies the relation 

J(-/3) = J(/3). 

In order that 0 be real the imaginary part of equation (A2) must be zero. Since 
the imaginary part of a complex number Z can be written as 9m Z = (Z - Z* )/2i, it 
follows from equation (A 2) that 0 will be real only if 



f D(0)J(/3)e i/3 ^ d/3 - f D* (0)J(0)e~^P^ d/3 = 0 

J -OO j -oO 


By changing the variable in the first integral from /3 to -/3 we obtain 
f D(-/3)J(-/3)e" W d(-/3) - f D* (/3) J(/3)e _i ^ d/3 = 0 

J 00 ■/ CO 


(A3) 


and since J(-/3) = J(/3), equation (A3) reduces to 


f ^ [d(-/ 3) - D*(/3)] J(/3)e“ i/3 ^ d/3 = 0 


Hence 0 will be real only if D(-/3) = D*(/3). 
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APPENDIX B 


C(k) 

C P 

D(a),D(j3) 

E(or),E(/3) 

F 

f 

H 

h 

9m 

a a 

M 

J ± 

K (<r,0) 

k 

k 

m 

k(or, j3) 

N 

P 

P 

Q 

^s 

‘Is 

Re 

s > s 0 

T 

t 

U 


SYMBOLS 

function defined in eq. (21) 
specific heat of fluid 

unknown in complex form of integral equation 

unknown in real form of integral equation 

function in specified temperature distribution (5) 

defined from F by eq. (32) 

function in specified pressure distribution (6) 

defined from H by eq. (31) 

integrals defined in eq. (39) 

imaginary part of a function 

unit vectors in X- and Y-directions 

functions defined in eq. (51) 

kernel defined by eq. (38) 

Fourier transform variable 
thermal conductivity of porous material 
real kernel defined by eq. (44) 
temperature parameter, (t^ - t^)/(t^ - t^) 
pressure parameter, (pg - Pj)/(Pq - Pj) 
pressure 

dimensionless heat flux, Q s 5A m (t^ - t^) 

conduction heat flux at surface 

real part of a function 

upper and lower surfaces of porous wall 

dimensionless temperature, (t - ^/(tj - t^) 

temperature 

dimensionless velocity, - u 

K (P Q " P X ) 



<1 


Darcy velocity 

dimensionless normal velocity at upper surface, ^ - v 

k (P 0 -Pi) S 

normal velocity leaving upper surface 

complex potential, W = ^ + i <p 

dimensionless coordinate, x/6 

dimensionless coordinate, y/6 

rectangular coordinates 

complex variable, X + iY 

transform variable 

transform variable 

defined in eq. (37) 

defined in eq. (43) 

wall thickness 

dependent variable defined in eq. (24) 
permeability of porous material 

P c n *(P 0 ~ Pi) 
parameter, — t — — 

2k n 

m r 

fluid viscosity 
fluid density 

parameter in pressure and temperature variations, eqs. (53) and (54) 
potential, imaginary part of W, (p Q - p)/(p Q - Pj) 
real part of W 
dimensionless gradient 


pressure variation 
along upper surface 
temperature variation 
at lower boundary of wall 

at large negative and positive x on upper surface 
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2 

OO 

Superscript: 

* 


value on upper surface having maximum deviation from value at large |x 
coolant reservoir conditions 


complex conjugate 
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